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Abstract 70-kDa Heat shock proteins are ATP-driven molecular chaperones that perform a 
myriad of essential cellular tasks. Although structural and biochemical studies have shed some light on 
their functional mechanism, the fundamental issue of the role of energy consumption, due to 
ATP- hydro lysis, has remained unaddressed. Here we establish a clear connection between the non- 
equilibrium nature of Hsp70, due to ATP hydrolysis, and the determining feature of its function, namely 
its high affinity for its substrates. Energy consumption can indeed decrease the dissociation constant 
of the chaperone-substrate complex by several orders of magnitude with respect to an equilibrium 
scenario. We find that the biochemical requirements for observing such ultra-affinity coincide with the 
physiological conditions in the cell. Our results rationalize several experimental observations and pave 
the way for further analysis of non-equilibrium effects underlying chaperone functions. 
DOI: 10.7554/el_ife.0221 8.001 
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Introduction 

ATP-driven molecular chaperones play a central role in protecting cells against proteins that could 
unfold or misfold because of mutations, various stresses or fluctuations and ultimately result in cyto- 
toxic aggregates (Bukau et aL, 2006; Hartl et aL, 2011). 70-kDa Heat shock proteins (Hsp70s) stand 
out for several reasons: they are possibly the most ubiquitous, they function as monomers, and they 
supervise a plethora of diverse cellular processes {Mayer and Bukau, 2005] Zuiderweg et aL, 2013) 
such as protein translation {Kramer et aL, 2009), protein trafficking (M attack et aL, 1999; Neupert 
and Brunner, 2002), the disassembly of protein complexes {Booking et aL, 2011), signaling {Pratt 
and Toft, 2003) and protein degradation {Hohfeld et aL, 2001). All these tasks crucially depend on 
the high-affinity binding of Hsp70s to substrate proteins during a complex ATP-driven conformational 
cycle {Mayer, 2013). The ATP- and ADP-bound states of Hsp70 (Hsp70-ATP and Hsp70-ADP respec- 
tively) and their interconversion play a major role in the chaperone functional cycle: the nature of the 
bound nucleotide affects the affinity of the chaperone for its substrates, with Hsp70-ADP binding 
the substrate more stably than Hsp70-ATP {Schmid et aL, 1994; Theyssen et aL, 1996; Packschies 
et aL, 1997; Gisler et aL, 1998; Russell et aL, 1998; Laufen et aL, 1999; Mayer et aL, 2000). 

Intriguingly, several experimental evidences suggested that the effective affinity of Hsp70 for sub- 
strates when the chaperone was running through its ATP-hydrolysis driven cycle was significantly 
higher than both the ones of Hsp70-ATP and Hsp70-ADP {Laufen et aL, 1999; Wittung-Stafshede 
et aL, 2003). Because this remarkable result would not be possible within the boundaries of thermo- 
dynamic equilibrium, it is therefore necessary to clarify how ATP hydrolysis, and thus energy consump- 
tion, can affect the binding strength of Hsp70s to their substrates. 

Results and discussion 

According to the consensus Hsp70 cycle {Figure 1) substrate binding/unbinding takes place with rates 
that depend on the state of the bound nucleotide {k™ TPl kf TPl k™ DPl kf DP ) {Schmid et aL, 1994; Gisler 
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eLife digest Proteins perform numerous essential tasks in cells. Most of these tasks require the 
protein to have a very specific structure, which is maintained by a balance of chemical and physical 
interactions. However, this delicate balance is vulnerable to excessive heat, changes in the pH of 
the cell, and certain chemicals. As a consequence, proteins could lose their specific structure and 
stop working. 

Cells employ a group of specialized proteins — called chaperones — to check that other proteins 
have the correct structure, and to 'refold' those that do not. The Hsp70 chaperone family needs 
energy to do its job, and it gets this energy from a molecule called ATP. However, the exact way 
that Hsp70s work and use this energy is not fully understood. 

One major puzzle is how Hsp70 binds to a protein to fold it up. Previous experiments suggested 
that this binding is particularly effective if Hsp70 can adopt different structures as part of a complex 
cycle governed by ATP. Now, De Los Rios and Barducci reveal that the energy released from breaking 
down ATP molecules enables this extra-efficient binding to occur. According to the proposed model, 
this is possible under some conditions that are normally found in cells. These include having many 
more Hsp70 proteins than target proteins, and producing energy at extremely high rates from ATP. 
The specific kinetic properties of the different structures Hsp70 can form are also crucial. 

More generally, the principle that energy consumption enhances binding could be 
extended beyond chaperone proteins and represent a general mechanism for other biomolecular 
systems. 

DOI: 10.7554/el_ife.0221 8.002 



et a/., 1998; Mayer et a/., 2000), with Hsp70-ATP exchanging the substrate two to three orders of 
magnitude faster than Hsp70-ADP The conversion from Hsp70-ADP to Hsp70-ATP occurs through a 
nucleotide exchange process which is described here at an effective level as a simple first-order reac- 
tion with rate k DTl or ko T in the presence of a bound substrate (Figure 1, and 'Materials and methods' 
for a full derivation). The conversion from Hsp70-ATP to Hsp70-ADP can occur by means of two dif- 
ferent processes: nucleotide exchange (dashed arrows in Figure 1, with rate kj x Dl and k e -^ s in the pres- 
ence of the substrate) and ATP hydrolysis (red arrows in Figure 1), whose rate depends on the absence 
or presence of a bound substrate {k h and kl respectively). The total conversion rate from Hsp70-ATP to 
Hsp70-ADP is thus kjo = kj x D + k^ (and analogous expressions in the presence of a substrate). In the cell, 
several cochaperones tune the exchange and hydrolysis rates: J-domain proteins (JDPs) enhance the 
rate of ATP hydrolysis, and nucleotide exchange factors (NEFs) catalyze nucleotide release (Youker 
and Brodsky, 2007] Kampinga and Craig, 2010). Within the present description, cochaperones are 
not taken into account explicitely. Rather, their action is captured as a modulation of the cycle time- 
scales. In particular, JDPs are known to bind the substrate and subsequently interact with Hsp70, 
enhancing ATP hydrolysis. Consequently here only the hydrolysis rate in the presence of the substrate, 
kl, is affected by the action of JDPs. 

A significant difference between hydrolysis and nucleotide exchange must be outlined here: 
ATP-hydrolysis, at variance with exchange, results into a net production of ADP and a loss of ATP In a 
cellular perspective, the ATP and ADP concentrations are kept fixed by energy-consuming chemo- 
stats. In vivo, ATP hydrolysis is therefore an intrinsically non-equilibrium process. 

Because we aim here at elucidating the relation between energy consumption and substrate 
affinity, we determine the effective dissociation constant of the system, K effl which provides a coarse- 
grained measure of how well Hsp70s can bind their substrates through their cycle. /C e ff is defined in the 
usual way as K e ^ = [S][Hsp70]/[Hsp70-S], where [Hsp70] is the total concentration of Hsp70 not bound 
to a substrate ([Hsp70] = [Hsp70-ATP]+[Hsp70-ADP]), [Hsp70-S] is the total concentration of substrate- 
bound chaperone ([Hsp70-S] = [Hsp70-ATP-S]+[Hsp70-ADP-S]) and [S] is the concentration of free 
substrate. 

In the absence of hydrolysis, no energy is consumed and all the reactions of the Hsp70 cycle are 
at equilibrium. In this scenario, where all the reactions are driven by thermal fluctuations, the detailed 
balance rule holds and each branch of the biochemical cycle is individually balanced ('Materials 
and methods'). In fact, in this case, the ratio between the forward and backward rates for each reaction 



De Los Rios and Barducci. eLife 2014;3:e02218. DOI: 10.7554/eLife.02218 



2 of 10 



eLIFE 



Research article 



X ^ ADP 



^ ATP 



► 



off 
ATP 



Biochemistry | Biophysics and structural biology 



A A 
V , 




«. — .7 A > 






i^ex 1 
*TD | 




I, I Lrex^ 1 
*DT *TD I 


kg 


*DT 


r 


> 


r *"*""*\ 1 


> 


r 



.Off 



ix(ADP) _ ADP ix (ATP) 

"I. on D 

ADP 



l, off 
i/ (ATP) _*ATP 

D "k on 

ATP 



Figure 1. Canonical Hsp70 biochemical cycle. The model takes into account four states in Hsp70 (NBD is schemati- 
cally represented here in green, SBD in orange), which are defined by substrate binding and by the nature of the 
bound nucleotide (ADP or ATP). The rates of the substrate binding/unbinding process (horizontal blue lines) are 
influenced by the nucleotide (k™ TP , k°^ TP , k™ TP , k°^ TP ). ADP-bound states are converted to ATP-bound states through a 
nucleotide exchange process (vertical solid blue lines) with rates k DT , k s DT . The ATP to ADP conversion can occur by 
means of either a nucleotide exchange process (dashed blue lines) with rates kj X D , k^f or ATP-hydrolysis (red lines) 
with rates k hf k*. 
DOI: 10.7554/eLife.0221 8.003 



is completely determined by the free energy difference between the two states, for example 
k s DT jkj D = exp (G Hsp70 . ADP . s - G H s P 7o-atp-s)- The equilibrium effective dissociation constant {K e J { ) can be 
easily determined as 

IS K (ATP) ,S K (ADP) 
j^eq _ _DT D TD D 



DT + TD 

where K { * TP) and K { * DP) are the dissociation constants of the Hsp70-ATP-substrate and Hsp70-ADP- 
substrate complexes, respectively. Not surprisingly, corresponds to a weighted average of K { £ TP) 
and K { * DP \ and it cannot be lower than K { £ DP) , usually the lowest of the two (Schmid et aL, 1994; Gisler 
et aL, 1998] Mayer et aL, 2000). As a matter of fact, in vivo would be close to its typical upper 
bound, namely K { ^ TP \ because of the cellular excess of ATP over ADP The equilibrium dissociation 
constant can be in principle measured in experiments where Hsp70 hydrolysis deficient mutants are 
used {McCarty and Walker, 1991), at varying ratios of the concentrations of the two nucleotides. 

At variance with equilibrium, when ATP hydrolysis is turned on, the energy budget along the cycle 
is not restricted to the free-energy differences between the different states. Rather, the dissipated 
energy E diss , which is related to ATP hydrolysis, must be taken into account. Thus detailed balance is 
broken and pairwise reactions cannot be solved individually as in the equilibrium case. 

Nonetheless, even in non-equilibrium conditions a steady-state solution of the cycle exists ('Materials 
and methods'), and it provides an expression for the non-equilibrium dissociation constant {K^). 

j^neq _ \ DT ATP ^ TD ADP ^ ATP ADP ) \ DT TD / \ K DT K ATP K ADP ^ TD ADP ATP J\_ \ 

eff k (k off +k s +/r s \k on + k (k off + k s 4- k s \k on +(k 5 + k s \k on k on 

DT \^ADP ^ DT ^ TD ) ATP ^ TD \ ATP ^ DT ^ TD ) ADP ^ \^DT ^ TD ) ATP ADP 1° \ 

In this scenario, ATP hydrolysis is controlled by the fixed basal hydrolsysis rate k h and the substrate- 
enhanced rate k^, which is further modulated in cellular conditions by JDP co-chaperones. The ratio 
k^/kb, that measures the hydrolysis acceleration, is thus a natural parameter to characterize the beha- 
vior of the system. 
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In order to prove the intimate relation between this quantity and the total energy consumption, we 
report in Figure 2A the experimentally-measurable hydrolysis flux, P diss , defined as 

Pdiss = K [Hsp70 ■ ATP] + k s h [Hsp70 ■ ATP • S] 

as a function of k^/k^. Here we consider the experimentally determined parameters for the Escherichia 
coli DnaK-DnaJ system (see Table 1), and concentrations that roughly mimic cellular conditions 
([Hsp70] tot = 40 [jM and substoichiometric substrate, here [S] tot = 4 pM). Not surprisingly, the hydrolysis 
flux increases with k^/k^, before saturation, corresponding to a regime dominated by the rate-limiting 
exchange process {Gassier et al., 2001). 

In Figure 2B we report K™ q , as a function of the acceleration ratio k^/k^. As k^/k^ increases and more 
energy is consumed, K™ q decreases, until it becomes lower than K { £ DP) by several orders of magnitude. 
Non-equilibrium conditions lead thus to a dramatic increase of the affinity of Hsp70s for their sub- 
strates, that could not be achieved at equilibrium, where the effective dissociation constant would be 
bounded between the ones of the ATP-bound and ADP-bound states. We dub such effect ultra-affinity 
in analogy with energy-consuming ultrasensitivity observed in many enzymatic systems {Goldbeter 
and Koshland, 1981). 

The lower bound of the non-equilibrium dissociation constant is /C^ q = k^ p /k™ TP , which can be 
achieved for extremely high values of k s h lkh (Figure 2B). This regime corresponds to the limiting case 
of substrate binding exclusively to Hsp70-ATP, which has the fastest binding rate, and being released 
exclusively from the ADP-bound state, which has the slowest unbinding rate. Our analysis indicates 
that this theoretical limit, recently hinted at {Zuiderweg et al., 2013), likely pertains to a regime that 
is not accessible to Hsp70s. 

It must be stressed here that this lower bound, as well as ultra-affinity, depends on the kinetic prop- 
erties of the cycle and does not rely on the dissociation constants of any nucleotide-bound state. In 
order to better elucidate this point, we move beyond the experimentally measured rates for the DnaK/ 
DnaJ system, and we explore the theoretical dependence of K™ q on the time-scale separation between 
the binding/unbinding kinetics in the two states. To this aim, in Figure 3A we report K™ q as a function 
of both k^/kb and the ratio /cadp/^atp/ which measures the time-scale separation in the unbinding 
kinetics, while keeping K { ^ DP) unchanged. For k^p ^ ^atp the non-equilibrium dissociation constant is 
bound to the equilibrium range. Ultra-affinity can be achieved only for /c A d P //c A t P <1, and is more pro- 
nounced for larger time-scale separations. In the limit /c A d P « k^ P , both the binding and the unbinding 
processes of the ADP-state become negligible and the non-equilibrium dissociation constant reduces to 

-(ATP) ^DT (^DT + ^TD + ) 

D k ik s + k ex,s + k s \ 

DT \ DT TD ^ K h ) 

which is reported in Figure 2B (dashed line) and provides a good approximation of the exact behavior 
in the physiologically-accessible range of k s h lkh. 

Strikingly, ultra-affinity depends also on the concentration of the substrate as can be inferred from 
the explicit expression for K™ q . This is another effect intrinsically tied to the non-equilibrium nature of 
the cycle. Indeed, at equilibrium, dissociation constants do not depend on the total species concen- 
tration since they simply encode the difference of free-energy between the bound and unbound 
states. We thus explore in Figure 3B the dependence of K™ q on both k^/k^ and the ratio between the 
total substrate and the total chaperone concentrations (for [Hsp70] tot = 40 fjM). Two distinct regimes 
can be observed here with a sharp transition occuring at [S] tot /[Hsp70] tot = 1 . In the excess of chaperone 
([S] to t/[Hsp70] tot <1) we observe ultra-affinity and the previously described behavior, whereas in the 
excess of substrate ([S] tot /[Hsp70] tot >1) the gain in affinity as a function of k^/k^ is limited and K™ q never 
exceeds its equilibrium range. This effect can be easily rationalized considering that in the latter con- 
dition, substrate binding to ADP-bound state becomes dominant and the system cannot exploit the 
time-scale separation in binding/unbinding kinetics to achieve ultra-affinity. 

Quite surprisingly, the experimental characterization of the Hsp70-substrate dissociation constant 
in physiological, non-equilibrium conditions is extremely limited. However, ultra-affinity was implicitely 
suggested in a series of works assessing the binding of different substrates to Hsp70s, always in the 
presence of a co-localized JDP, thus ensuring maximal hydrolysis acceleration upon substrate binding 
(Misselwitz et al., 1998, Laufen et al., 1999, Sullivan et al., 2001). In all these assays the substrate 
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Figure 2. Effect of ATP-hydrolysis on K^ q . Total energy consumption (A) and effective non-equilibrium dissociation 
constant of the Hsp70-substrate complex (B) is plotted as a function of the hydrolysis acceleration ratio k s h /k h , for 
the DnaK/DnaJ/substrate system with concentrations [Hsp70] tot = 40 uM and [S] tot = 4 uM (see 'Materials and 
methods' for the parameters), The approximate dissociation constant K^ q 0 is also plotted for comparison (black 
dashed line). The green region comprised between K { ^ TP) and K { ^ DP) corresponds to the range of affinities accessible 
at equilibrium (no hydrolysis). The red-to-yellow region corresponds to the values of the dissociation constants that 
are exclusively accessible to the non-equilibrium regime. The region where red fades to yellow (10 3 < k s h /k h < 10 4 ) 
corresponds to the transition from physiological to non-physiological values of hydrolysis acceleration. 
DOI: 10.7554/el_ife.0221 8.004 



was observed to bind more efficiently in the presence of ATP, that drives the chaperone through its 
cycle, than in the presence of ADP, that instead blocks the system in the Hsp70-ADP state, the one with 
the lowest thermodynamic dissociation constant. We know of only one case where ultra-affinity has 
been carefully measured for the Hsp70 system interacting with a substrate peptide fused to a J-domain 
(Wittung-Stafshede et a/., 2003). The dissociation constant observed with this setup («0.22 nM) was 
two orders of magnitude smaller than both the measured K [ £ TP) and K [ £ DP) (>30 nM). If we combine 
our model with the specific rates provided in {Wittung-Stafshede et al., 2003) and values of k s h lkh 
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Table 1. Parameters of the model compatible to what reported in the literature, we 

obtain a predicted in the range 0.15-0.6 nM, 
which is in excellent agreement with the experi- 
mental value. All these findings strongly suggest 
that ultra-affinity becomes manifest when sub- 
strate binding is coupled with enhanced hydrol- 
ysis acceleration by the colocalization of the 
Hsp70 binding region and of a J-domain. 

Our analysis of the cycle has unveiled the 
conditions that Hsp70 must satisfy to exhibit 
ultra-affinity: (/) the substrate-exchange rates of 
the ADP-bound state must be significantly slower 
than the ATP-state so that extremely different 
timescales can be exploited; (/'/) ATP hydrolysis 
must be enhanced by orders of magnitude to 
fully enter the ultra-affinity regime; (///) the chaper- 
one must be in excess over the substrate. All these 
conditions are typically met in the cell by canon- 
ical Hsp70 chaperones, such as bacterial DnaK 
and cytosolic Hsc70 in eukaryotes: Hsp70-ADP is 
known to have limited exchange kinetics; 
Hsp70s are known to work only in partnership 
with JDPs; Hsp70s are highly abundant and typi- 
cally in excess over JDPs {Finka and Goloubinoff, 2013). Cellular conditions seem thus to be optimal 
for ultra-affinity. 

Our results provide an additional example of how ATP hydrolysis can be exploited by cells to over- 
come the constraints set by equilibrium thermodynamics. The key role of energy consumption in driving 
biochemical cycles for performing a variety of functions is well established (Schnakenberg 1976; Hill, 
2005; Ge et al, 2012) and recently it has been recognized in cellular processes such as sensing, sign- 
aling and adaptation (Q/an and Reluga, 2005, Lan et al, 2012, Mehta and Schwab, 2012). 

Notably, the ultra-affinity concept proposed here shares some similarities with the well-known 
kinetic proofreading {Hopfield, 1974) since in both cases chemical energy consumption is used to 
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Parameters used in the model, from various sources. 
Notable cases are: 

*This corresponds to the reaction Hsp70ATP-DnaJ 2 -S -»• 

Hsp70-ADP-DnaJ 2 -S. 

DOI: 10.7554/el_ife.0221 8.005 




Figure 3. Dependence of K^ q on time-scale separation and on stoichiometric ratio. (A) Non-equilibrium dissociation constant as a function of the 
hydrolysis acceleration ratio k s h /k h and of the time-scale separation between the ATP- and ADP-state, expressed as the ratio between the substrate 
unbinding rates between the ATP- and ADP-state. (B) Non-equilibrium dissociation constant as a function of the hydrolysis acceleration ratio lc^/lc h and 
of the stoichiometric ratio between the total substrate and Hsp70 concentrations. The color codes are the same as in Figure 2, green for the region 
accessible in equilibrium, and red-to-yellow for the region accessible in non-equilibrium. The blue line is the non-equilibrium dissociation constant 
reported in Figure 2B. 
DOI: 10.7554/el_ife.0221 8.006 
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increase the binding affinity of specific molecules beyond their equilibrium value. However, a signifi- 
cant dissimilarity can be outlined: kinetic proofreading exploits multiple non-equilibrium steps to 
enhance a pre-existing difference in binding affinity among various substrates, whereas ultra-affinity is 
achieved due to the specific ability of the substrate to induce upon binding a non-equilibrium transi- 
tion in the receptor {i.e., the chaperone). 

We expect that a similar approach may be applicable to more specialized Hsp70s, such as HscA 
(Hesterkamp and Bukau, 1998), and to other chaperones, with possibly more complex cycles, such 
as the GroEL/S system, HsplOOs and Hsp90s. As well, it would not be surprising to discover other 
molecular machines working along the same principles, so that the present study might also provide 
a general scheme applicable beyond Hsp70 chaperones. 

Materials and methods 

We frame here the model in terms of deterministic mass-balance equations, which describe the evolu- 
tion of the concentrations of the different states. This treatment is appropriate for the Hsp70 system 
since both chaperones and substrates are highly abundant in standard cellular conditions. In the case 
of less abundant molecular species fluctuations might not be negligible. Therefore a more detailed, 
stochastic description of the process would be required (Ge et al., 2012). 

Coarse graining of the nucleotide exchange processes 

The exchange process corresponds to the reactions 

k-ATP ^ADP 

Hsp70 • ATP + ADP Hsp70 + ATP + ADP Hsp70 • ADP + ATP 

<— <- 

^ATP k ADp 

The corresponding rate equations are 

d[Hs P 70-ATP] = _ k _^ [Hsp7Q ATp] + ^ [ATp][Hsp70] 

d[H J 7Q] = " {k + ATP [ATP] + k + ADP [ADP]) [Hsp70] + k ATP [Hsp70 ■ ATP] 
+ k- ADP [Hsp70-ADP] 

d[Hs P 70-ADP] = _^ [Hsp7Q ADp] + ^ [ADp] [Hsp?o] 

where k ATP and k ADP are the ATP and ADP release rates, respectively, and k ATP and k ADP are the ATP and 
ADP binding rates, respectively. 

In steady state we can obtain the concentration of the nucleotide-free state from the second equa- 
tion, and substitute it in the equations for the nucleotide-bound state concentrations, obtaining, for 
example (from the first equation) 

d [ H5p7 ,°- ATP ] = -k- TP - r ^ DP ] THspyO-ATP] 
oft ATP /c; TP [ATP] + /c; DP [ADP] L J 



+ k- ADP k ATP [ATP] -[Hsp70-ADP] 

ADP /c; TP [ATP] + /c; DP [ADP] L J 



and analogously for the equation for [Hsp70-ADP]. 

As a consequence, the expressions for the effective exchange rates kj^ and k^ T are 

/.ex _ /c^ Dp [A DP] 
K TD ~ K ATP~ 



fcA7p[ATP] + k* [ADP] 
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k + ATP [AJP] 



K DJ K ADP 



^ T p[ATP] + ^ DP [ADP] 



Typically the nucleotide release rates are much slower than the nucleotide binding rates {Table 1), and 
the above expressions correctly capture that nucleotide release sets upper bounds for the rates, which are 
then further modulated by the partitioning between ATP and ADP binding. Importantly, the equilibrium 
constant between the ATP-bound and ADP-bound states is correctly reproduced by the ratio kj^/k^j. 
Analogous expressions can be obtained for the exchange process in the presence of a bound substrate. 

Coarse graining of the co-chaperone action 

The acceleration of hydrolysis by JDPs is described here as a modulation of k s h in the range k^k^kf, ^, 
where /ch max is the maximal experimentally determined JDP/substrate accelerated hydrolysis rate. In this 
work we extend our analysis also to larger values of k s h for completeness. A more complete description of 
the underlying co-chaperone binding/unbinding reactions would only overburden the present model with- 
out providing further insights. The dimerization of DnaJ into DnaJ 2 has also been considered as implicit. 

Parameters of the model 

The rates that we have used to solve the model equations and for the data in Figure 2,B have been 
taken from studies of the DnaK/DnaJ/GrpE system, consistently with previous modeling (Hu et al., 
2006), and are reported in Table 1. The value of k JD ex,s has been obtained from the other parameters, 
using the relation 

K (ATP) ,ex,S 
i ex,S _ i ex fX D DT 
TD ~ TD ^(ADP) ^ex 

which holds for a cycle at thermodynamic equilibrium because the free energy difference accumulated 
over a cycle is AG = 0 (Ge et al., 2012). The upper bound for the hydrolysis acceleration induced by 
JDPs {kl max) corresponds to the hydrolysis rate experimentally observed in saturation of DnaJ concom- 
itant with the presence of the substrate {Laufen et al., 1999). 

The rate of synthesis of ATP from ADP (the microscopic reverse of hydrolysis) is neglected here 
because it was experimentally proven to be below the level of detectability even in the presence of 
the hydrolysis-accelerating JDP cochaperones (Russell et al., 1998), that is /c synth <1 0~ 6 s _1 . 

In all the calculations the chaperone concentration is 40 pM. The ratio [ATP]/[ADP] = 10 has been 
used throughout the calculations, approximately matching the physiological ratio. 

Equilibrium solution of the Hsp70 cycle 

At equilibrium the detailed balance rule holds, which implies that the scheme in Figure 1 can be 
solved by balancing each branch of the cycle individually. K^J corresponds to the solution of the sys- 
tem of equations 

kZp [S][Hsp70 • ATP] = k° A % [Hsp70 • ATP • S] 



k e TD [Hsp70 • ATP] = k DT [Hsp70 • ADP] 
k° A n DP [S] [Hsp70 • ADP] = k° A ff DP [Hsp70 • ADP • S] 
k 5 h [Hsp70 • ATP • S] = k 5 DT [Hsp70 • ADP • S] 

[S] = [S] tot - ([Hsp70 • ATP • S ]+[ Hsp70 • ADP • S]) 

where [S] tot is the total substrate concentration. 

Non-equilibrium solution of the Hsp70 cycle 

In non-equilibrium the detailed balance rule does not hold anymore, and individual branches do not 
lead to a solution of the full cycle. Yet, the steady-state solution still exists and can be found solving 
the steady-state mass-balance equations for the scheme in Figure 1, namely 
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kZ P [S] [Hsp70 • ATP] - (k s h + kf TP ) [Hsp70 • ATP • S] + k s DT [Hsp70 -ADP-S] = 0 

kf TP s [Hsp70 • ATP • S] - (k h + /c° n TP [S]) [Hsp70 • ATP] + k DT [Hsp70 • ADP] = 0 
k °AD P [S] [Hsp70 - ADP] - (k s DT + kf DP ) [Hsp70 • ADP • S] + k 5 h [Hsp70 • ATP • S] = 0 
kfj [Hsp70 - ADP-S] - (k DT + k° A n DP [S]) [Hsp70 • ADP] + k h [Hsp70 • ATP] = 0 
[S] = [S] tot - ([Hsp70-ATP-S]+ [Hsp70- ADP-S]) 
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